Towards clinical application: prompt gamma imaging of passively scattered proton fields with a knife-edge slit camera.

Prompt γ-ray imaging with a knife-edge shaped slit camera provides the possibility of verifying proton beam range in tumor therapy. Dedicated experiments regarding the characterization of the camera system have been performed previously. Now, we aim at implementing the prototype into clinical application of monitoring patient treatments. Focused on this goal of translation into clinical operation, we systematically addressed remaining challenges and questions. We developed a robust energy calibration routine and corresponding quality assurance protocols. Furthermore, with dedicated experiments, we determined the positioning precision of the system to 1.1 mm (2σ). For the first time, we demonstrated the application of the slit camera, which was intentionally developed for pencil beam scanning, to double scattered proton beams. Systematic experiments with increasing complexity were performed. It was possible to visualize proton range shifts of 2-5 mm with the camera system in phantom experiments in passive scattered fields. Moreover, prompt γ-ray profiles for single iso-energy layers were acquired by synchronizing time resolved measurements to the rotation of the range modulator wheel of the treatment system. Thus, a mapping of the acquired profiles to different anatomical regions along the beam path is feasible and additional information on the source of potential range shifts can be obtained. With the work presented here, we show that an application of the slit camera in clinical treatments is possible and of potential benefit.

Analytical computation of prompt gamma ray emission and detection for proton range verification.

A prompt gamma (PG) slit camera prototype recently demonstrated that Bragg Peak position in a clinical proton scanned beam could be measured with 1-2 mm accuracy by comparing an expected PG detection profile to a measured one. The computation of the expected PG detection profile in the context of a clinical framework is challenging but must be solved before clinical implementation. Obviously, Monte Carlo methods (MC) can simulate the expected PG profile but at prohibitively long calculation times. We implemented a much faster method that is based on analytical processing of precomputed MC data that would allow practical evaluation of this range monitoring approach in clinical conditions. Reference PG emission profiles were generated with MC simulations (PENH) in targets consisting of either (12)C, (14)N, (16)O, (31)P or (40)Ca, with 10% of (1)H. In a given geometry, the local PG emission can then be derived by adding the contribution of each element, according to the local energy of the proton obtained by continuous slowing down approximation and the local composition. The actual incident spot size is taken into account using an optical model fitted to measurements and by super sampling the spot with several rays (up to 113). PG transport in the patient/camera geometries and the detector response are modelled by convolving the PG production profile with a transfer function. The latter is interpolated from a database of transfer functions fitted to MC data (PENELOPE) generated for a photon source in a cylindrical phantom with various radiuses and a camera placed at various positions. As a benchmark, the analytical model was compared to MC and experiments in homogeneous and heterogeneous phantoms. Comparisons with MC were also performed in a thoracic CT. For all cases, the analytical model reproduced the prediction of the position of the Bragg peak computed with MC within 1 mm for the camera in nominal configuration. When compared to measurements, the shape of the profiles was well reproduced and agreement for the estimation of the position of the Bragg peak was within 2.7 mm on average (1.4 mm standard deviation). On a non-optimized MATLAB code, computation time with the analytical model is between 0.3 to 10 s depending on the number of rays simulated per spot. The analytical model can be further used to determine which spots are the best candidates to evaluate the range in clinical conditions and eventually correct for over- and under-shoots depending on the acquired PG profiles.

Dose ratio proton radiography using the proximal side of the Bragg peak.

PURPOSE: In recent years, there has been a movement toward single-detector proton radiography, due to its potential ease of implementation within the clinical environment. One such single-detector technique is the dose ratio method in which the dose maps from two pristine Bragg peaks are recorded beyond the patient. To date, this has only been investigated on the distal side of the lower energy Bragg peak, due to the sharp falloff. The authors investigate the limits and applicability of the dose ratio method on the proximal side of the lower energy Bragg peak, which has the potential to allow a much wider range of water-equivalent thicknesses (WET) to be imaged. Comparisons are made with the use of the distal side of the Bragg peak. METHODS: Using the analytical approximation for the Bragg peak, the authors generated theoretical dose ratio curves for a range of energy pairs, and then determined how an uncertainty in the dose ratio would translate to a spread in the WET estimate. By defining this spread as the accuracy one could achieve in the WET estimate, the authors were able to generate lookup graphs of the range on the proximal side of the Bragg peak that one could reliably use. These were dependent on the energy pair, noise level in the dose ratio image and the required accuracy in the WET. Using these lookup graphs, the authors investigated the applicability of the technique for a range of patient treatment sites. The authors validated the theoretical approach with experimental measurements using a complementary metal oxide semiconductor active pixel sensor (CMOS APS), by imaging a small sapphire sphere in a high energy proton beam. RESULTS: Provided the noise level in the dose ratio image was 1% or less, a larger spread of WETs could be imaged using the proximal side of the Bragg peak (max 5.31 cm) compared to the distal side (max 2.42 cm). In simulation, it was found that, for a pediatric brain, it is possible to use the technique to image a region with a square field equivalent size of 7.6 cm(2), for a required accuracy in the WET of 3 mm and a 1% noise level in the dose ratio image. The technique showed limited applicability for other patient sites. The CMOS APS demonstrated a good accuracy, with a root-mean-square-error of 1.6 mm WET. The noise in the measured images was found to be σ = 1.2% (standard deviation) and theoretical predictions with a 1.96σ noise level showed good agreement with the measured errors. CONCLUSIONS: After validating the theoretical approach with measurements, the authors have shown that the use of the proximal side of the Bragg peak when performing dose ratio imaging is feasible, and allows for a wider dynamic range than when using the distal side. The dynamic range available increases as the demand on the accuracy of the WET decreases. The technique can only be applied to clinical sites with small maximum WETs such as for pediatric brains.

Absolute prompt-gamma yield measurements for ion beam therapy monitoring.

Prompt-gamma emission detection is a promising technique for hadrontherapy monitoring purposes. In this regard, obtaining prompt-gamma yields that can be used to develop monitoring systems based on this principle is of utmost importance since any camera design must cope with the available signal. Herein, a comprehensive study of the data from ten single-slit experiments is presented, five consisting in the irradiation of either PMMA or water targets with lower and higher energy carbon ions, and another five experiments using PMMA targets and proton beams. Analysis techniques such as background subtraction methods, geometrical normalization, and systematic uncertainty estimation were applied to the data in order to obtain absolute prompt-gamma yields in units of prompt-gamma counts per incident ion, unit of field of view, and unit of solid angle. At the entrance of a PMMA target, where the contribution of secondary nuclear reactions is negligible, prompt-gamma counts per incident ion, per millimetre and per steradian equal to (124 ± 0.7stat ± 30sys) × 10(-6) for 95 MeV u(-1) carbon ions, (79 ± 2stat ± 23sys) × 10(-6) for 310 MeV u(-1) carbon ions, and (16 ± 0.07stat ± 1sys) × 10(-6) for 160 MeV protons were found for prompt gammas with energies higher than 1 MeV. This shows a factor 5 between the yields of two different ions species with the same range in water (160 MeV protons and 310 MeV u(-1) carbon ions). The target composition was also found to influence the prompt-gamma yield since, for 300/310 MeV u(-1) carbon ions, a 42% greater yield ((112 ± 1stat ± 22sys) × 10(-6) counts ion(-1) mm(-1) sr(-1)) was obtained with a water target compared to a PMMA one.

Prompt gamma imaging of proton pencil beams at clinical dose rate.

In this work, we present experimental results of a prompt gamma camera for real-time proton beam range verification. The detection system features a pixelated Cerium doped lutetium based scintillation crystal, coupled to Silicon PhotoMultiplier arrays, read out by dedicated electronics. The prompt gamma camera uses a knife-edge slit collimator to produce a 1D projection of the beam path in the target on the scintillation detector. We designed the detector to provide high counting statistics and high photo-detection efficiency for prompt gamma rays of several MeV. The slit design favours the counting statistics and could be advantageous in terms of simplicity, reduced cost and limited footprint. We present the description of the realized gamma camera, as well as the results of the characterization of the camera itself in terms of imaging performance. We also present the results of experiments in which a polymethyl methacrylate phantom was irradiated with proton pencil beams in a proton therapy center. A tungsten slit collimator was used and prompt gamma rays were acquired in the 3-6 MeV energy range. The acquisitions were performed with the beam operated at 100 MeV, 160 MeV and 230 MeV, with beam currents at the nozzle exit of several nA. Measured prompt gamma profiles are consistent with the simulations and we reached a precision (2σ) in shift retrieval of 4 mm with 0.5 × 10(8), 1.4 × 10(8) and 3.4 × 10(8) protons at 100, 160 and 230 MeV, respectively. We conclude that the acquisition of prompt gamma profiles for in vivo range verification of proton beam with the developed gamma camera and a slit collimator is feasible in clinical conditions. The compact design of the camera allows its integration in a proton therapy treatment room and further studies will be undertaken to validate the use of this detection system during treatment of real patients.

Real-time proton beam range monitoring by means of prompt-gamma detection with a collimated camera.

Prompt-gamma profile was measured at WPE-Essen using 160 MeV protons impinging a movable PMMA target. A single collimated detector was used with time-of-flight (TOF) to reduce the background due to neutrons. The target entrance rise and the Bragg peak falloff retrieval precision was determined as a function of incident proton number by a fitting procedure using independent data sets. Assuming improved sensitivity of this camera design by using a greater number of detectors, retrieval precisions of 1 to 2 mm (rms) are expected for a clinical pencil beam. TOF improves the contrast-to-noise ratio and the performance of the method significantly.

Prompt gamma imaging with a slit camera for real-time range control in proton therapy.

Treatments delivered by proton therapy are affected by uncertainties on the range of the beam within the patient, requiring medical physicists to add safety margins on the penetration depth of the beam. To reduce these margins and deliver safer treatments, different projects are currently investigating real-time range control by imaging prompt gammas emitted along the proton tracks in the patient. This study reports on the feasibility, development and test of a new concept of prompt gamma camera using a slit collimator to obtain a one-dimensional projection of the beam path on a scintillation detector. This concept was optimized, using the Monte Carlo code MCNPX version 2.5.0, to select high energy photons correlated with the beam range and detect them with both high statistics and sufficient spatial resolution. To validate the Monte Carlo model, spectrometry measurements of secondary particles emitted by a PMMA target during proton irradiation at 160 MeV were realized. An excellent agreement with the simulations was observed when using subtraction methods to isolate the gammas in direct incidence. A first prototype slit camera using the HiCam gamma detector was consequently prepared and tested successfully at 100 and 160 MeV beam energies. Results confirmed the potential of this concept for real-time range monitoring with millimetre accuracy in pencil beam scanning mode for typical clinical conditions. If we neglect electronic dead times and rejection of detected events, the current solution with its collimator at 15 cm from the beam axis can achieve a 1-2 mm standard deviation on range estimation in a homogeneous PMMA target for numbers of protons that correspond to doses in water at the Bragg peak as low as 15 cGy at 100 MeV and 25 cGy at 160 MeV assuming pencil beams with a Gaussian profile of 5 mm sigma at target entrance.

Water equivalent path length measurement in proton radiotherapy using time resolved diode dosimetry.

PURPOSE: To verify water equivalent path length (WEPL) before treatment in proton radiotherapy using time resolved in vivo diode dosimetry. METHODS: Using a passively scattered range modulated proton beam, the output of a diode driving a fast current-to-voltage amplifier is recorded at a number of depths in a water tank. At each depth, a burst of overlapping single proton pulses is observed. The rms duration of the burst is computed and the resulting data set is fitted with a cubic polynomial. RESULTS: When the diode is subsequently set to an arbitrary depth and the polynomial is used as a calibration curve, the "unknown" depth is determined within 0.3 mm rms. CONCLUSIONS: A diode or a diode array, placed (for instance) in the rectum in conjunction with a rectal balloon, can potentially determine the WEPL at that point, just prior to treatment, with submillimeter accuracy, allowing the beam energy to be adjusted. The associated unwanted dose is about 0.2% of a typical single fraction treatment dose.